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We derive new constraints set by SNIa experiments ('gold' data sample of Riess et at), X-ray 
galaxy cluster data (Allen et al. Chandra measurements of the X-ray gas mass fraction in 26 clus- 
ters), large scale structure (Sloan Digital Sky Survey spectrum) and cosmic microwave background 
(WMAP) on the quartessence Chaplygin model. We consider both adiabatic perturbations and in- 
trinsic non-adiabatic perturbations such that the effective sound speed vanishes (Silent Chaplygin). 
We show that for the adiabatic case, only models with equation of state parameter \cc\ < 10~ 2 are 
allowed: this means that the allowed models are very close to ACDM. In the Silent case, however, 
the results are consistent with observations in a much broader range, —0.3 < a < 0.7. 



I. INTRODUCTION 

Two of the major puzzles of contemporary cosmology 
are the nature of dark matter and of dark energy, the two 
big players in the cosmic arena. So far, the only knowl- 
edge of these components refers to their density fraction 
and to their equation of state, and even on these numbers 
we still have a large uncertainty. Greater still is the igno- 
rance of their clustering properties: although we assume 
that essentially all the dark matter clusters in observable 
objects and all the dark energy does remain quite ho- 
mogeneous, this is to a large extent only a simplification 
rather a consequence of observations. 

It is therefore no surprise that many works are cur- 
rently devoted to the possibility of merging the two puz- 
zles into a single one: that is, finding a single origin for 
both dark matter and dark energy. A possibility is to 
assume an interaction between dark matter and dark en- 
ergy or to fit both into a single complex field Q. 
However, these models still contain two separate fields 
that account ultimately for the two components. On a 
different level lies the hypothesis that there is a single 
fluid that behaves as dark matter or dark energy accord- 
ing to the background or the local density. Since there 
is only one unifying dark-matter-energy component, be- 
sides baryons, photons and neutrinos, this model is usu- 
ally referred to as quartessence 0|. A phenomenologi- 
cal prototype of quartessence models is the generalized 
Chaplygin model 0. 0. E3|. an exotic fluid with an in- 
verse power law homogeneous equation of state (EOS), 
p = — Af 4 ( Q+1 )/p Q , where M has dimension of mass and 
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a > —1 is a dimensionless parameter (a = 1 is the orig- 
inal Chaplygin gas). Such equation of state leads to a 
component which behaves as dust in the past and as cos- 
mological constant in the future. For a — the model 
reduces to ACDM 0. 

For a wide range of values of the parameter a, the 
quartessence Chaplygin model is compatible with several 
cosmological tests that involve only the background met- 
ric Q. Nevertheless, problems may occur when one con- 
siders perturbations. For instance, the CMB anisotropies 
spectrum is strongly suppressed with respect to ACDM 
[2g. Further, it was shown that, unless a is very much 
close to a = 0, the mass pow er sp ectrum presents strong 
instabilities and oscillations pi]]. A quantitative CMB 
analysis 0] (see also 0] for a pre- WMAP analysis of the 
generalized Chaplygin gas as dark energy) found that the 
parameter a should be small (a < 0.2 at 95% CL), even 
including an additional CDM component, and smaller 
still without the latter. We will show below that in fact 
this limit reduces to \a\ < 0.01. It is clear that this 
strong limit is due to the finite sound speed of the Chap- 
lygin fluid: when it is sufficiently large and positive, it 
prevents clustering on small scales and therefore intro- 
duces a cut in the power spectrum which is at odds with 
observations; when it is negative, on the other hand, the 
clustering at small scales is enhanced beyond control. 

It is important to stress, however, that unlike the back- 
ground tests, the perturbation analysis need further as- 
sumptions beyond the equation of state. In 0], it was 
shown that it is possible to avoid the mass power spec- 
trum problem if, for instance, pressure perturbation van- 
ishes. This can be done by introducing a special type 
of intrinsic entropy perturbation such that the effective 
sound speed of the cosmic fluid vanishes (we de- 
note this as silent perturbations) . We call the Chaplygin 
quartessence model with vanishing pressure perturbation 
Silent Chaplygin and refer to the standard case as adia- 
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batic Chaplygin. 

The main goal of this paper is to show that, unlike 
adiabatic Chaplygin, Silent Chaplygin is consistent with 
CMB data for a wide range of parameters. Besides CMB, 
we also consider consistency of Silent Chaplygin with cur- 
rent data from large scale structure, from type la super- 
novae (SNIa) and baryon fraction in galaxy clusters. Us- 
ing the latest data sets, we review all these tests using 
a careful treatment, presenting the outcome of a com- 
bined analysis of the data for both Silent and adiabatic 
Chaplygin. 



II, BASIC EQUATIONS 

In this section we give a brief and basic description of 
the model. The conservation equation for the generalized 
Chaplygin gas component in a Robertson-Walker metric 
is solved by 



Pch = Pchol-wo + 



1 



w 



,3(a+l) 



(1) 



1 ChO 



where a is the scale factor (a = 1 today) and w — 1 ' 

is the present equation of state. 

The equation-of-state parameter (w = p/p) and the 
adiabatic sound velocity (c 2 = dp/dp) for the Chaplygin 
gas, are given by 



wch{a) 



WQ 



-w + (1 + w )a~ 3(Q+1) 



and 



c 2 s Ch(a) = -awch(a). 



(2) 



(3) 



From the equation above, it is clear that at early times, 
when a — > 0, we have wch —> 0, and the fluid behaves as 
non relativistic matter. At late times, when a ^> 1, we 
obtain wch — > — 1- Further, the adiabatic sound speed 
has maximum value when the equation of state param- 
eter is minimum. Therefore, the epoch of cosmic accel- 
eration, in these models, coincides with that of large c 2 . 
As remarked in this is a characteristic of all unified 
models in which the equation of state is convex, i.e., is 

such that = -^=f < 0. However, this property is not 
mandatory for a general fluid. Models with concavity 
changing equations of state, may have c 2 negligibly small 
when the energy density reaches its minimum value. This 
is an important property that must be considered in con- 
structing acceptable adiabatic quartessence models (l^ . 

The perturbation equations for a fluid with equation 
of state w and adiabatic sound speed c 2 in synchronous 
gauge are (for the sake of simplicity, here we neglect 
baryons and radiation, and we assume that both spa- 
tial curvature and the anisotropic perturbation vanish) 

5' = -3{wT + c 2 s S - wS) - {1 + w){9 + ^) (4) 
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V = 



3(1 + w)n 2 on 

~^2k 2 
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where derivatives are with respect to logo. Here 5 is the 
density contrast, f2 is the density parameter, hh and r\ are 
metric perturbations, 9 is the divergence of the velocity 
perturbation, T is the entropy perturbation and 7i = 
da/dr, t being the conformal time. To these equations, 
we must add the equations for baryons and relativistic 
fluid. The assumption of a silent universe requires that 
131 



Sp = pT + c 2 s {Sp) = 



so that 



w 

We adopt therefore these equations with w = wch and 
Cs = CsCh- With this choice, T — aS and models with 
a = are both, silent and adiabatic. 



III. OBSERVATIONAL CONSTRAINTS 

In the subsections below, we will derive constraints on 
the parameters a and Wq from four data sets: SNIa, X- 
ray cluster gas fraction, SDSS power spectrum and CMB. 

A. SNIa 

The luminosity distance of a light source is defined in 
such a way as to generalize to an expanding and curved 
space the inverse-square law of brightness valid in a static 
Euclidean space, 



d L 



\AttT ) 



1/2 



(8) 



In Eq. JSl, L is the absolute luminosity and T is the 
measured flux. 

For a source of absolute magnitude M, the apparent 
bolometric magnitude m{z) can be expressed as 

m(z) =M-51ogifo+25 + 51ogZ>L 

=M + 42.3841 - blogh +5logD L , (9) 
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where Dl is the luminosity distance in units of H 1 (we 
are using c = 1), 



D L = f ' [E(z')]- l dz' , 
Jo 



(10) 



where 
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E z (z) = [u b h- 2 (1 + z) 3 
+ (1 - cj b h- 2 ) \-wo + (1 + w ) (1 + z) 3(1+Q) 



(11) 



Here, z is the redshift, uj b = f^ft 2 , and ft. is the Hubble 
constant in units of 100km/ sMpc . 

For the SNIa analysis, we use the 'gold' data set of 
Riess et al. The data in this sample is given in 

terms of the extinction corrected distance modulus, [Xq = 
rn — M. To obtain the likelihood of the parameters, we 
use a x 2 



,2 statistics such that, 



X {a,w ,h,LO b ) = 

157 



E 



jUo,i - 42.3841 + 5 log ft - 5 log /^(z^a, Wo, h,Ub) 



(12) 



are the estimated errors in the individual 
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distance moduli, including uncertainties in galaxy red- 
shift and also taking into account the dispersion in su- 
pernovae redshift due to peculiar velocities (see 0] for 
details) . To determine the likelihood of the parameters a 
and wo, we marginalize the likelihood function over h and 
tu b . We use two different priors in this work. We consider 
flat priors when combining the SNIa and clusters results 
with CMB, while in this and in the next subsection we 
adopt Gaussian priors such that ft = 0.72 ± 0.08 and 
Lo b = 0.0214±0.002 [19]. The results of our SNIa analysis 
are displayed in Fig. ^ In the figure, we show 68% and 
95% confidence level contours in the (a, wq) plane. 



B. X-ray Cluster Gas Fraction 

Clusters of galaxies are the most recent large-scale 
structures formed and the largest gravitationally bound 
systems known. Therefore, the determination of their 
matter contents is important since cluster properties 
should approach those of the Universe as a whole. By 
measuring the baryon mass fraction Sl b /Q m in rich clus- 
ters, and combining this ratio with fl b determinations 
from primordial nucleosynthesis, constraints on fl m can 
be placed [20]. Further, by assuming that baryon mass 
fraction in clusters of galaxies is independent of redshift, 
it is also possible to constrain the geometry and, con- 
sequently, the dark energy density. A method based on 
this idea was suggested by Sasaki [5]J and Pen j23| and 
further developed and applied by Allen, Schmidt, and 
Fabian El. 



SNIa 

h=0.72 ± 0.08 

co,=0.021 4 ±0.002 
a 




Figure 1: Constant confidence contours (68% and 95%) in the 
(a, wo) plane allowed by SNIa data, as described in the text 



In this section we use the new data set of Allen et 
al. Q to constraint the Chaplygin models. These 
authors extracted from Chandra observations the x-ray 
gas mass fraction (/ gas ) of twenty six massive, dynami- 
cally relaxed galaxy clusters, with redshifts in the range 
0.08 < z < 0.9, and that have converging / gas within 
a radius r25oo (radius encompassing a region with mean 
mass density 2500 times the critical density of the Uni- 
verse at the cluster redshift). 

To determine the confidence region of the parameters 
of the model, we use the following x 2 function in our 
computation: 



X 2 (a, w ,h,uj b ,b) 



26 

E 

i=l 



[fgas d ( z t> a , w o,h,uj b ,b) - f, 



gas,' 



(13) 

where z,, / gas ,ij and er/ gaB i are, respectively, the redshifts, 
the SCDM (ft = 0.5) best-fitting values, and the symmet- 
ric root-mean-square errors for the 26 clusters as g iven 
in 0. In Eq. /™° d is the model function 
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3/2 



(l + 0.19Vft)^ 



0.5 c£- 



(14) 



Here, d,A = (1 + z) _2 c?l is the angular diameter distance 
to the cluster, b is a bias factor that takes into account the 
fact that the baryon fraction in clusters could be lower 
than for the Universe as a whole and 



eff 



(1 - u b h- 2 ) (1 + Wo) 



u b h 



(15) 



is the effective matter density parameter 01 . In our 
computations, we first marginalize analytically over the 
bias factor assuming that it is Gaussian-distributed with 
b = 0.824 ±0.089 [21123. To determine the likelihood of 
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Figure 2: Constant confidence contours (68% and 95%) in the 
(a, wo) plane from cluster f gaa data as described in the text. 
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Figure 4: CMB likelihood functions for a and wo in the adi- 
abatic case. The horizontal lines mark the 68% and 95% c.L, 
top to bottom. 
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Figure 5: Constant confidence contours (68% and 95%) in the 
(a, wo) plane allowed by CMB, as described in the text 
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Figure 3: Constant confidence contours (68% and 95%) in the 
(a, Wo) plane allowed by CMB, as described in the text. 



the parameters a and wq, we next marginalize the like- 
lihood function over h and u>b- As remarked before, we 
assume here a Gaussian prior such that h = 0.72 ± 0.08 
and uj b = 0.0214 ± 0.002 0. 
In Fig. 12 we show the 68% and 95% confidence con- 
tours on the parameters a and wo determined from the 
Chandra data. 



C. CMB 

Here we compare the model to the combined tem- 
perature and polarization power spectrum estimated by 
WMAP [2f|. To derive the likelihood we adopt a ver- 



sion of the routine described in Verde et al. [27:], which 
takes into account all the relevant experimental proper- 
ties (calibration, beam uncertainties, window functions, 
etc). 

Our theoretical model depends on two Chaplygin pa- 
rameters, four cosmological parameters and the overall 
normalization N: 

a,Wo,n s ,h,LOb, r, N. (16) 

The overall normalization has been integrated out numer- 
ically. We calculate the theoretical Ci t t spectra by a mod- 
ified parallelized CMBFAST jjj3 code that includes the 
full set of perturbation equations [3, [[lj with the addi- 
tion of non- adiabatic pressure perturbations. We do not 
include gravitational waves and the other parameters are 
set as follows: T cmb = 2.726X, Y He = 0.24, N v = 3.04. 

We evaluated the likelihood on a unequally spaced grid 
of roughly 50,000 models (for each normalization) with 



5 




-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 
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Figure 6: CMB likelihood functions for a and tuo in the silent 
case. The horizontal lines mark the 68% and 95% C.I., top to 
bottom. 



the following top-hat broad priors: w E (—1,-0.5), 
a € (-0.8,1.5), n s € (0.8, 1.2), u>b G (0.02, 0.025), 
t € (0.,0.3). For the Hubble constant we adopted the 
top-hat prior h <E (0.6,0.8); we also employed the HST 
result f.18] h = 0.72 ± 0.08 (Gaussian prior). 

In Fig. 03 we show the confidence region on a, wq in the 
adiabatic case, after marginalization over all the other 
parameters with flat priors: as it can be seen, this is 
the most stringent test among those studied in this pa- 
per. As anticipated, it restricts a to be very close to 0: 
in other words, observations of CMB demands that the 
background of adiabatic Chaplygin be almost indistin- 
guishable from ACDM. In Fig. 0we contrast this result 
with the silent case: now the allowed region for a widens 
a lot, encompassing values close to unity. What is partic- 
ularly interesting is that even with a = 1, Silent Chap- 
lygin is consistent at the 2a level, while it was obviously 
ruled out in the adiabatic case. In Figs. 31 and we 
show the one-dimensional likelihood for a and wq, with 
marginalization over all the other parameters. 



D. Large scale structure 

The results of Sandvik et al. (2004) [HI have shown 
that the matter power spectrum of adiabatic Chaplygin 
is plagued by strong instabilities and oscillations for any 
a significantly different from zero, leading to a stringent 
upper limit to \a\. In it has been shown that these 
problems do not occur for the silent model. To test this in 
a more quantitative way we perform here a simplified like- 
lihood analysis with the following procedure. We com- 
pared the baryon spectra of the silent case to the mat- 
ter power spectrum of Sloan Digital Sky Survey (SDSS) 
as obtained in Tegmark et al. (2004) (2a], cutting at 
k = 0.02/i/Mpc, using the likelihood routine provided by 
M. Tegmark [23, and marginalizing over the amplitude 
(i.e., we are comparing only the slope of the spectrum, 
not its absolute value). In order to save computer time 
we restricted ourselves to a subset of parameters, fixing 
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Figure 7: SDSS likelihood functions at 68% and 95 % c.l. for 
a and wo in the silent case. 



LOb = 0.023, n s = 1 and h = 0.7. As we show in Fig. 
the results are rather similar to the CMB case. Since 
we used only a very small subset of the whole param- 
eters space we will not include the SDSS likelihood in 
the combined likelihood of the next section. The loss in 
constraining power is acceptable since one can easily see 
that the CMB and the SDSS confidence regions are over- 
lapping to a good extent. Moreover, extending the range 
of parameters as in the CMB case the SDSS constraints 
would become weaker and would not add much to the 
information from the CMB. 

We also calculated the baryon fluctuation variance as 
(i.e. the absolute normalization of the spectrum) for the 
same subset of parameters. In Fig. [HI we show the con- 
tour plots of the likelihood assuming a Gaussian distri- 
bution of the observed as with as = 0.9 ± 0.1. This test 
is meant to be only qualitative since the value of as that 
is derived from data is strongly degenerate with Sl m and 
with the bias factor and / or it depends on calibration ob- 
tained with A-body simulations based on some specific 
cosmological model. We only notice that for the Silent 
Chaplygin case there exists a large region in our param- 
eter space which accounts for values of as which are con- 
sistent with the standard estimation. It is nevertheless 
interesting to note that the as analysis shows preference 
for lower values of wq (for a fixed a) than the CMB test. 
This could help to further reduce the parameters space. 



IV. COMBINED CONSTRAINTS AND 
CONCLUSION 

Here we present a combined analysis of the constraints 
discussed in the previous section. In Fig. El we display 
the allowed region of the parameters wq and a from a 
combination of data from SNIa, clusters and CMB in 
adiabatic Chaplygin. The constraints on the parameter 
a are roughly the same as those obtained from CMB 
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Figure 8: Likelihood function for as at 68% and 95 % c.l. for 
a and wo in the silent case. 
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Figure 9: Constant confidence contours (68% and 95%) in 
the (a, wq) plane allowed by SNIa, clusters and CMB, as de- 
scribed in the text. 



alone. Only models with a very close to the ACDM limit 
are allowed. Although including SNIa and clusters in 
the CMB analysis almost do not affect the constraints 
on a, they tighten the constraints on the parameter Wq, 
reducing even more the allowed region of the parameters 
space in the adiabatic case. The final result (95% c.l.) 
for the adiabatic case is 



-0.005 < a < .01, -0.8 < w Q < -0.7. 



(17) 



In contrast, in Fig. the Silent Chaplygin model is 
consistent with the observables considered here for a wide 
range of parameters. The constraints in the silent case 
are (95% c.l.) 
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CMB+SNIa+Clusters 



a 



Figure 10: Constant confidence contours (68% and 95%) in 
the (a, wo) plane allowed by SNIa, clusters and CMB, as de- 
scribed in the text 
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Figure 11: Likelihood functions for a and wo in the adiabatic 
case for SNIa, clusters and CMB. The horizontal lines mark 
the 68% and 95% c.l., top to bottom. 
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(18) 



Figure 12: Likelihood functions for a and wo in the silent 
case for SNIa, clusters and CMB. The horizontal lines mark 
the 68% and 95% c.l., top to bottom. For comparison, we 
plot as dotted lines the combined likelihood functions of the 
adiabatic case. 
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with a remarkable fifty-fold extension in the range of a. 
The one-dimensional likelihood functions are shown in 

Figs. El and C3 

The idea of unifying dark matter and dark energy 
through a single component has motivated many works 
in the last few years. Most of these investigations 
concentrated their efforts in analyzing the generalized 
Chaplygin fluid, that is considered the protot ype of the 
quartessence models. After the results of jfl ll(lllll |. it be- 
came clear that adiabatic Chaplygin, as quartessence, is 
ruled out unless the parameter a is very close to zero. In 
the present paper we confirmed this result by using cur- 
rent CMB data, extending the quantitative results of ^lj 
to quartessence models. In fact, for an effective sound 
speed different from zero, any quartessence model with 
a convex EOS will suffer the same kind of problem 
There are two simple ways to get around this problem: 
choosing an EOS not convex [16] or introducing entropic 
perturbations as we did in this paper. This shows that 
indeed quartesence models of the Chaplygin type (one 
fluid, two parameters) can be considered as real alterna- 



tives to ACDM (two fluids, one parameter). It is clear 
however that these phenomenological models should be 
further investigated: the challenge is to connect them to 
a more fundamental theory. 

Finally, we remark that one possible source of problems 
for Silent Chaplygin is lensing skewness 31]. It would 
be remarkable if higher-order or non-linear effects would 
prove necessary to rule out, or strongly constrain, these 
models. 



Acknowledgments 

We thank Maurlcio Calvao, Roberto Colistete, Martin 
Makler and Ribamar Reis for useful discussions. The 
CMB computations have been performed at CINECA 
(Italy) under the agreement INAF@CINECA. We thank 
the staff for support. IW is partially supported by the 
Brazilian research agency CNPq. 



[1] L. Amendola, Phys. Rev. D 62, 043511, (2000). 
[2] R. Mainini & S.A. Bonometto, Phys. Rev. Lett. 93, 
121301 (2004). 

[3] M. Makler, S.Q. Oliveira & I. Waga, Phys. Lett. B 555, 
1 (2003). 

[4] A. Kamenshchik, U. Moschella & V. Pasquier, V., 2001, 

Phys. Lett. B 511, 265 (2001). 
[5] N. Bilic, G.B. Tupper & R.D. Viollier, Phys. Lett. B 535, 

17 (2002). 

[6] M.C. Bento, O. Bertolami & A. A. Sen, Phys. Rev. D 66 

, 043507 (2002). 
[7] P.P. Avelino, L.M.G. Beca, J. P.M. de Carvalho & 

C.J.A.P. Martins, JCAP 09, 002 (2003). 
[8] M. Makler, S.Q. Oliveira & I. Waga, Phys. Rev. D 64, 

123521 (2003); A. Dev, D. Jain & J.S. Alcaniz, Astron. 

Astrophys. 417, 847 (2004); R. Colistete Jr. & J.C. Fab- 

ris, Class. Quant. Grav. 22, 2813 (2005); M.C. Bento, O. 

Bertolami, N.M.C. Santos, A. A. Sen, Phys. Rev. D 71, 

063501 (2005). 

[9] D. Carturan & F. Finelli, Phys. Rev. D 68, 103501 
(2003). 

[10] H. Sandvik, M. Tegmark, M. Zaldarriaga & I. Waga, 

Phys. Rev. D 69, 123524 (2004). 
[11] L. Amendola, F. Finelli, C. Burigana & D. Carturan, 

JCAP 07, 005 (2003). 
[12] R. Bean & O. Dore, Phys. Rev. D 68, 023515 (2003). 
[13] R.R.R. Reis, I. Waga, M.O. Calvao & S.E. Joras, Phys. 

Rev. D 68, 061302(R) (2003). 
[14] W. Hu, Astrophys. J. 506, 485 (1998). 
[15] R.R.R. Reis, M. Makler & I. Waga, Class. Quant. Grav. 



22, 353 (2005). 

[16] M. Makler, R. R. R. Reis, L. Amendola & I. Waga, in 
preparation 

[17] A.G. Riess et al, Astrophys. J. 607, 665 (2004). 
[18] W. Freedman et al, Astrophys. J. 553, 47 (2001). 
[19] D. Kirkman et al, Astrophys. J. Suppl., 149, 1 (2003). 
[20] S.D.M. White & C. S. Frenk, Astrophys. J. 379, 52 

(1991); S.D.M. White, J.F. Navarro, A.E. Evrard & C. 

Frenk, Nature (London) 366, 429 (1993). 
[21] S. Sasaki, Publ. Astron. Soc. Jpn., 48, L119 (1996). 
[22] U. Pen, New Astronomy 2, 309 (1997). 
[23] S.W. Allen, R.W. Schmidt & A.C. Fabian, Mont. Not. 

R. Astron. Soc. 334, Lll (2002). 
[24] S.W. Allen et al, Mont. Not. R. Astron. Soc. 353, 457 

(2004). 

[25] D. Rappeti, S.W. Allen & J. Weller, Mont. Not. R. As- 
tron. Soc. 360, 555 (2005). 

[26] G. Hinshaw et al. [WMAP collaboration], Astrophys. J. 
Suppl. 148, 135 (2003). 

[27] L. Verde et al. [WMAP collaboration], Astrophys. J. 
Suppl. 148, 195 (2003). 

[28] M. Tegmark, et al. [the SDSS collaboration], Phys. Rev. 
D 69, 103501 (2004); M. Tegmark, et al. [the SDSS col- 
laboration], Astrophys. J. 606, 702 (2004). 

[29] http:/ /space. mit.edu/home/tegmark/sdss. html 

[30] U. Seljak & M. Zaldarriaga, Ap. J. 469, 437 (1996). 

[31] R.R.R. Reis, M. Makler & I. Waga, Phys. Rev. D 69, 
101301(R) (2004). 



